Integrable open supersymmetric U model with boundary impurity 



A. Foerster 1 *, K.E. Hibberd 2 t, J.R. Links 3 * and I. Roditi 4 . § 
(February 1, 2008) 

^nstituto de Fisica da UFRGS, 
Av. Bento Goncalves 9500, Porto Alegre - RS, Brazil. 

o" 

, 2 ' 4 Centro Brasileiro de Pesquisas Fisicas, 

Rua Dr. Xavier Sigaud 150, 22290-180, Rio de Janeiro - RJ, Brazil. 

3 Department of Mathematics, 



3 

o 



c3 



in 



The University of Queensland, 4072, Australia. 



An integrable version of the supersymmetric U model with open boundary conditions and an 
impurity situated at one end of the chain is introduced. The model is solved through the nested 
algebraic Bethe ansatz method so that the Bethe ansatz equations are obtained. 



PACS numbers: 71.20.Fd, 75.10.Jm, 75.10.Lp 



I. INTRODUCTION 

T3 

The study of one-dimensional models of correlated electrons in the context of Luttinger liquids attracts a wealth of 
activity presently. Particularly, such models in the presence of defects or impurities are of particular interest in the 
i i ' effort to unify the results of both theoretical and experimental investigations. As a result, substantial research has 
| been devoted to this field from a variety of approaches. 

Lee and Toner first considered the problem of an impurity coupled to a Luttinger liquid Jl[] and the Kondo (magnetic 
impurity) problem has been studied in detail using conformal field theory J^-Q], perturbation theory j^], renormal- 
ization group H and Bethe ansatz methods [f-10|. Still the problem of strongly correlated electron models with 



' impurities remains not well understood and recent discrepancies between observed experimental results and theoreti- 
, cal predictions Jll]] have justified the search for a deeper understanding. 

Some of the well studied one-dimensional integrable spin models with impurities include the (supersymmetric) t-J 
0^ [ model [|l2]-[l5| and Heisenberg chain [f7] JTl| , p^ -[l9| . However one discovers that when impurities are introduced into a 
periodic lattice the Hamiltonian, though integrabl e, may contain terms that are physically unrealistic as the system 
suffers from an absence of backscattering ^|,|ll|-^3 15 ll| . 



There exists another class of one-dimensional models with impurities. These are systems constructed with open 
boundary conditions where the impurity appears on the end of the chain. It is desirable to introduce impurities into 
the system in such a way that integrability is maintained. The quantum inverse scattering method (QISM) allows for 



such a construction and following Sklyanin's approach |20 one may obtain an integrable impurity model with open 
boundary conditions. This has been demonstrated for a number of models with Kondo impurities 21 -|2^]. 

This article presents a new boundary impurity for the supersymmetric U model p8[ maintaining integrability. Since 
its discovery, the supersymmetric U model has received considerable attention due partly to the fact that the model 
• »"H , contains a free integrability preserving coupling parameter. The Bethe ansatz equations and thermodynamic aspects 
of the model on the closed chain were discussed in (29) and an anisotropic generalization was presented in |jo| . 

The boundary impurity we propose for the supersymmetric U model, is constructed from an i?-matrix solution 
of the Yang-Baxter relation. The solution is obtained from the tensor product of one parameter family of typical 
four-dimensional irreducible representations of gl(2\l) characterized by distinct parameters. This produces boundary 
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terms with a different coupling of the Hubbard interaction, compared with that of the bulk Hamiltonian. In addition, 
the boundary term also contains a spin interaction which is absent in the bulk interactions. 

The paper is organised as follows. The construction of the Hamiltonian is presented in the next section. In section 
III the nested algebraic Bethe ansatz solution is presented and concluding remarks follow in section IV. 



II. THE SUPERSYMMETRIC U MODEL WITH BOUNDARY IMPURITY 



In this section, we introduce the basic ingredients to construct an integrable impurity model with open boundary 
conditions through the QISM. We begin by writing the i?-matrix |3l| 

R^iu) = ^^-p - p 2 + (iLi) 

— u — a — p u — a — p — 2 

where the a, arise from the one-parameter family of four-dimensional representations of g/(2| 1) and the projectors 
Pi, P2 and P 3 are given explicitly in the appendix. In particular, for (3 = a we recover the i?-matrix used to derive 
the usual supersymmetric U model p8j. 

This i?-matrix acts in the tensor product of two four-dimensional modules V a ® Vp and is a solution to the Z2-graded 
quantum Yang-Baxter equation 

R^(m ~ u 2 )R^{ui - u 3 )R%g(u 2 - u 3 ) = R%£(u 2 ~ u 3 )R^{m - u 3 )R^{u x - u 2 ). {11.2) 

The usual notation is adopted, that is, Rjk(u) acts on j-th and k-th superspaces and as an identity on the remaining 
superspace. Throughout the article we shall use the graded-tensor product law defined by 

(a ®i))(c®d) = ( ac ® bd). {11.3) 

This i?-matrix also satisfies unitarity 

Ri 2 (u)R 2 i(-u) = p(u), 

and crossing unitarity 

R%{~u + l)R£{u)=p{u), 

where p(u), p(u) are some scalar functions. 

In order to deal with open boundary conditions we follow the necessary supersymmetric generalization of Sklyanin's 
approach p0| and define the boundary transfer matrix t{u) as 

t(u) = str (^ + (u)T(u) J ft:_(u)T- 1 (-u)) . (II.4) 

Here T{u) is the monodromy matrix for a L-site chain 

T{u) = R$?(u)Rg?(u) ■ ■ ■ R^^R^iu), (II.5) 

K- (u) , K+ {u) both obey the reflection equations 

Rii{ui - u 2 )K-{ui)R 2 i(ui + u 2 )K^{u 2 ) = K-{u 2 )Ri 2 {ui + u 2 )K-(ui)R 2 i(ui - u 2 ), 
R^i 1 lst2 {-ui + u 2 )Kf {ui)Ri 2 {-ui -u 2 + V )K^ ( U2 ) 

= KX t2 {u 2 )R 2 i{-ui -u 2 + r 1 )K^{ui)R£ ^(- U i+u 2 ). {11.6) 

Above 77 is the so-called crossing parameter and sti stands for the supertransposition taken in the i-th space, whereas 
isti is the inverse operation of sU. Here we will choose K_(u), K + {u) to be both equal to the identity matrix. Note 
that in the expression for the monodromy matrix T(u) (II. 5) the presence of an i?-matrix intertwining between two 
different four-dimensional representations of <?Z(2|1) characterized by a and (3. The corresponding transfer matrix will 
be associated to a supersymmetric U m odel with an impurity in the boundary. 

Using the Yang-Baxter algebra (H.2), together with the reflection equations ( |H.6| ), it can be shown that the model 
is integrable, that is 

[r(ui),T(« 2 )] = 0. (II.7) 
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From this transfer matrix one may write the Hamiltonian for an open supersymmetric U model with a boundary 
impurity through the second derivative of t(u) with respect to the spectral parameter u, at u = 0. This is the 
procedure adopted when the supertrace of K + (u) is null [B2| . The Hamiltonian is given by 



H 



where 



q a+l _ g -q-l 

i H ' 

Ing 

r"(0) 
4(Ai + 2Si) 

h ~ X 1 i 

i=l 
+ 



2(A 1 + 2B 1 ) 



str ( K+ (O)Glo ) +2 str ( k' + (0)H* o ) + str ( K+ (0)(ff£ ) 2 



(II. 



A x = str X;(0), B x = str ( K+ (0)tf£ o , , 



In the above, -Pj,i+i denotes the graded permutation operator acting on qu antum spaces i and + With our choice 
for the boundary K matrices equal to the identity, the Hamiltonian fllLq ) simplifies to 



JT* = 



r 



(0) 



L-l 



Wi^y = g H ^ + i^W) [stro(GLo) + stro (( ^ o)2)] • (IL9) 

In view of the grading, the basis vectors of the module V can be identified with the eletronic states as follows 

|1) = | U) = cjcjlO) , |2) = | T) = c||0) , |3) = I 1) = c||0) , |4) = |0). (11.10) 

We adopt the notation that (cjj.) are spin up annihilation (spin down creation) operators. Also let rii^ = c\ a Ci ir7 
denote the number operator for electrons with spin a on site % and let the total number of electrons be rii — n^f + 
The spin operators S , £t , S z , 

= 4,T C U ' 5 'l = 4i c ^T: = o^.i _n i.T): 

form an sl(2) algebra. These operators commute with the Hamiltonian of the supersymmetric U model which is given 
by 

#i,i+l = - Y] {c\ c i+l a + h.c.) ( ) + - [n iA n lA + n i+1A n l+ld ] 

»=u V « / 

H — (cj cf Cj+i.-ffCj+i^ + h.c.) + (m + n i+1 ). 
a 

The impurity Hamiltonian constructed from (BT9) is found to be 



L-l 



H=^H i>i+1 +Hl , (11.11) 



where the boundary term is given 



h l,o = + X]( c i-i^ CL ' cr + k- c -)[^ + / 9 )"L-i,-o-("-L,- C r - 1) + a)(n L _i._ cr - l)n L _ a 

(7 

+mnL-i,-anL-a] - |(c^_ li<T cJ / _ lj _ .Ci ) _ er ci, )<7 + /i.e.) - /(S^^Sl + h.c.) + f ^ n L -i, a nL,-<T 

a 

-(a + + p)n L -i + (1 + a)n L ] - 4[/3(l + f3)n L ^n L -i,i + a ( x + a ) n L^n L ,i}}, (H.12) 
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with 



m = v/^(2 + a + /3) 2 - {a + 0) V(l + a)(l + /3), d(a,j3) = -y/T+~P[4y/Z + {a + /3) 2 (V^ - VTTa)], 
f={a-/3) 2 , e = 4^/3(1 + a)(l + /3), t = 4(a + /?) 2 (2 + a + /?), 

-409+1) 



(a + /5) 2 V(l + «)(l+/3), Z 



(a + /3) 2 (2 + a + /3) 2 



At this point a few interesting characteristics of the model sho uld be pointed out. In contrast to the usual super- 
symmetric U model, one may observe that the above Hamiltonian (11.11) contains an extra spin interaction term from 
the boundary contribution. Also the free parameter a appears on the boundary as well as in the bulk terms of the 
Hamiltonian, in contrast the Kondo impurity model presented in p7f . Furthermore, in the limit j3 — ► a we recover 
the usual supersymmetric U model p8[. 



III. THE BETHE ANSATZ SOLUTIONS 



In order to find the Bethe ansatz equations of the Hamiltonian (11.11), we use Sklyanin's generalized Bcthe ansatz 
method to treat open boundary conditions together with Babujian and Tsvelick's j34j approach for higher spin chains. 
Therefore, we introduce the following doubled auxiliary monodromy matrix 

U{u) = (f (u)k-(u)f-\-u)) , (HI.l) 

where we will chose K^{u) to be the identity matrix and T(u) is the monodromy matrix for a L-site chain 

T («) = R^{u)R^{u) ■ ■■Rg {L _ 1) (u)R% L (u). 
Above the matrix R a (it) acts on W ® V a , where W is a three-dimensional module and it is given by |3(| 

\ 



(HI.2) 
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(III.3) 



The doubled auxiliary monodromy matrix (III.l) acts on W ® V a 

U{u) = 



L-l 



(IH.4) 

Vp and can be represented by the matrix 



^ii(u) A 12 (u) B x {u) 
A 2 i{u) A 2 2{u) B 2 (u) 
d{u) C a (u) V{u) 



It is possible to show that it fulfills a modified Yang-Baxter equation of the form 

Ri 2 (ui - u 2 )U(ui) R 2 \(ui + u 2 )U{u 2 ) = U(u 2 ) R\ 2 (ui + u 2 )U(ui) R 2 \(ui - u 2 ), 



(HI.5) 



(HI.6) 
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where R(u) is the usual t — J i?-matrix acting on W ® W p5|,p6[ 
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(ni.7) 



From the doubled auxiliary monodromy matrix (111.1) we define an auxiliary transfer matrix 

f(u) = str (k + (u)U(u)) , (111.8) 

where K + (u) is chosen to be the diagonal matrix with ele ment s { — 1, —1, 1}. 

It can be shown that this auxiliary transfer matrix f (u) ( HI.S| ) commutes with the transfer matrix t(u) (II. 4), which 
means that they have a common set of eigenvectors. These vectors will be determined by applying the algebraic Bcthc 
ansatz to t(u). For this purpose we begin by writing the Bcthe vector |fi) (according to the first-level Bethe ansatz) 
as 

|0) =C ll (u 1 )---C lN (u N )\y)F^- t ». (111.9) 

is the pseudovacuum and the coefficients F ll '" lN will be determined later by the second level Bethe ansatz. The 
action of the elements of U(u) (III.5) on | W) is given by 

„ i fu+a+2^ Lt 



(u + /3 + 2) 



B d {u)\V) 

Att(«)|*) 
Add(u)\^) 



= 0, 



u + a + 2 

C d (u)\V) f 0, 

L 



(-u + P + 2) 



l*>, 



1 



— u + a + 1 
-u + a + 2 

—U + a + 1 
-u + a + 2 

u + a + 2 



L-l 



l*>- (b^d) 



L-l 



2u + 1 V — u + a - 



L-l 



(-U + /3 + 1) 
(-u + (3 + 2) 

(-U + /3 + 1) 
(-U + /3 + 2) \2u+ 1 

(u + (3 + 2) 



2u 



(-u + 13 + 2) 



(111.10) 



Using the Yang-Baxter relation (|III.6[ ) together with the expressions for the matrices wwR( u ) ( |IH.7| ) and U(u) 
(III. 5), we can obtain the commutation relations between T>(u), Adb(u) and Cd(u). These are complicated relations 
but they will simplify considerably if instead of using Adb(u) we work with Adb(u) according to the transformation 

A b d(u) = A hd (u) - -^—5 b dV(u). (III.ll) 
2u + 1 

This transformation will allow us to more easily recognize the "wanted" and "unwanted" terms. More explicitly we 
have 



A~bd(ui)C c (u 2 ) = 



T^(ui)Cb(u 2 ) 



(ui - u 2 - l)(ui + u 2 ) eb , , ihl , n ( ,r, , 

7 rr^T TTT r gh (ui + u 2 + 1 rJ3(«i - u 2 ) C e (u 2 )A gi (u 1 ) 

(Ul - u 2 )(ui + u 2 + 1) y 

"7 rr^ 2 — 77 7 r 9b A2 Ul + 1) CJ Ul )V(u 2 ) 

(«i + U2 + l)(2«i + l)(2«2 + l) cd ' 9V ' V ' 

+ 7 -7775 — 77 TT rf d (2 Ul + 1) C g ( Ul )A ic (u 2 ), 

(ui - u 2 )(2ui + 1) 

(ui - u 2 - l)(ui + u 2 ) „ , s t 2u 2 



(ui - u 2 )(ui + u 2 + 1) 

— — —C d (ui)Adb(u 2 ) 

ui + u 2 + 1 



C b (u 2 )V( Ul ) + 



(ui - u 2 )(2u 2 + I) 



Cb(ui)T>(u 2 ) 



(111.12) 



(111.13) 
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Here the matrix r(it), which satisfies the Yang-Baxter equation, takes the form, 

i v 

r b b l(u) = l, r h b j{u) = -, r b d t(u) = -, (b ± d, b,d= 1,2). (111.14) 

u — 1 u — 1 



By applying the transfer matrix f (it) ( HI. 8 ) to the Bethe vector |f2) ( III.9| ) we can solve the eigenvalue problem 

f(u)|fi) = A(«)|n). 

We find the eigenvalue to be 

Ar )- 2u ~ 1 f u + a + 2 \ L1 f u + l3 + 2 \ TT + - 1) 

9?/ 1 \ — it 4- /v 4- 9 / \ —7/ 4-/3 4-9 / 11 



2u + 1 V -u + a + 2 J \-u + P + 2 J £1 ( u _ Uj .)( u + Uj - + l) 

^ (u + % -)(,- % -lh (1) 

A ^ (it — UjjtU + Uj + 1) 



2u ( u+\\ L ( -it + a + \\ L 1 {-U + /3 + V S 



2it + l\ ~u J \-u + a + 2 J \-u + j3 + 2) ( u - Uj ){ 
provided the parameters {uj} satisfy a first Bethe ansatz equation given by 



—Uj \ L ( iij + a + 2 \ L 1 / Uj + (3 + 2 



ij + 1 J \ —Uj + a + 1 J \ —Uj + (3 + 1 _ 

"TT i u J -Vm + l){u 3 +V m + |) + 2ltj + 1 py (Uj - Itfc ) (Uj + lt fc + 1) pj (ttj - D m ~ \){Uj + V m - |) ^ 

Above, A^^it; {itj}) is the eigenvalue of the transfer matrix f W (it) for the reduced problem which arises out of the 
r(it) matrices from the first term in the right hand side of (III. 12). The reduced transfer matrix fW(it) may be 
recognized as that of the TV-site inhomogeneous XXX spin-i open chain and may be diagonalized following Ref. |2Gfj 
with the reduced boundary K matrices are equal to the identity. We find the eigenvalue of the transfer matrix (it) 
to be given by 

r ,s _ 2u-l t^t (n - t> m + f )(u + v m + |) 

2U + 1 pr (it - Uj)(u + Uj +1) yr (it - U m - |)(lt + U m - f ) 

2it 1 = 1 (it - itj - l)(zt + itj) 11^ (it - v m + ±)(it + v m - \) ' 
provided the parameters {v m } satisfy a second Bethe ansatz equation given by 

pr (fm - Uj - f )(n m + Uj - |) = pj- (i' m - v k - l)(v m + v k - 2) , m 

J — 1 k—1 

k y£ m. 



IV. CONCLUSION 

We have presented the supersymmetric U model with a new type of boundary impurity. The model is constructed 
with open boundary conditions and formulated within a modified QISM so that integrability is not lost. The impurity 
is introduced through an i?-matrix solution of the Yang-Baxter equation, built from the tensor product of two different 
f our-d imensional representations of gl(2\l) which are characterized by different parameters. The resulting Hamiltonian 



(11.11) contains an extra spin interaction term that does not appear in the usual supersymmetric U model. Another 
feature of this impurity model is that the free parameter of the bulk appears in the boundary terms, in contrast to 
the Kondo impurity model presented in [ ^7| . In the limit f3 —> a we recover the supersymmetric U model [ p8| on an 
open chain. In this work, we have also derived the Bethe ansatz equations using the nested algebraic Bethe ansatz 
method, which was appropriately modified to deal with a graded underlying algebra and open boundary conditions. 
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APPENDIX: PROJECTION OPERATORS FOR THE ^"-MATRIX 



We present the projectors for the tensor product V a ® Vp. For details on the construction of i?-matrices from these 
projectors, the reader is referred to |plj|. 



The tensor product basis is constructed from (H.1C) and the tensor product decomposition is V a ® Vp = V(Q, 0, \a + 
0) © V{— 1, —1, \a + (3 + 2) V(0, — 1, \a + + 1). It is not necessary to calculate the three projectors as one may use 
the relation I dent = ^2iPi- The simplest projectors to calculate are P\ and P3 which are given as follows 
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d x = l/[a(a + 1) + 2ap + f3(f3 + 1)], d 2 = l/(a + 

fa = V«/3(a + 1)08 + 1), fi = aVP(a + l), f 2 = + 



(IV.l) 
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